Langevin equations for competitive growth models 
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Langevin equations for several competitive growth models in one dimension are derived. For 
models with crossover from random deposition (RD) to some correlated deposition (CD) dynamics, 
with small probability p of CD, the surface tension f and the nonlinear coefficient A of the associated 
equations have hnear dependence on p due solely to this random choice. However, they also depend 
on the regularized step functions present in the analytical representations of the CD, whose expansion 
coefficients scale with p according to the divergence of local height differences when p — >■ 0. The 
superposition of those scaling factors gives u ^ p^ for random deposition with surface relaxation 
(RDSR) as the CD, and ~ p, A ~ p^^^ for ballistic deposition (BD) as the CD, in agreement 
with simulation and other scaling approaches. For bidisperse ballistic deposition (BBD), the same 
scaling of RD-BD model is found. The Langevin equation for the model with competing RDSR and 
BD, with probability p for the latter, is also constructed. It shows linear p-dependence of A, while 
the quadratic dependence observed in previous simulations is explained by an additional crossover 
before the asymptotic regime. The results highlight the relevance of scaling of the coefficients of step 
function expansions in systems with steep surfaces, which is responsible for noninteger exponents 
in some p-dependent stochastic equations, and the importance of the physical correspondence of 
aggregation rules and equation coefficients. 

PACS numbers: 81.15Aa, 05.40.-a, 05.50.+q, 68.55.-a 



I. INTRODUCTION 

The technological applications of thin films and multi- 
layers motivated intense theoretical study of growth mod- 
els in the last decades^"— . Many processes show evidence 
of a competition between different aggregation dynamics. 
For instance, this occurs in deposition of diamond-like 
carbon by plasma, where growth is mainly due to aggre- 
gation of slow radicals but ion bombardment is essential 
to create sp^ bonds'*. When physico-chemical conditions 
are continuously changed, such as in cyclical electrodepo- 
sition/dissolution of metals, competing dynamics are also 
present^. Consequently, many competitive growth mod- 
els were already proposed, with microscopic aggregation 
rules representing the atomistic dynamics. They are usu- 
ally defined on lattices, such as those with aggregation of 
different species of particles^"— and those mixing different 
aggregation rules for the same species^i"— . They usually 
show crossover effects from one dynamics at small times 
t and short length scales L to another dynamics at long 
t or large L, and in special cases anomalous roughening 
is presentSSii^. 

A widely studied group of models is that showing 
crossover from random (uncorrelated) deposition (RD) to 
some correlated deposition (CD) process, hereafter called 
RD-CD models. For RD-CD in general, aggregation of 
each incident particle follows the rules of the CD process 
with probability p and those of RD with probability 1—p 
(other models might also show the same crossover with a 
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parameter that is not a probabilitj^) . Another group of 
relevant models show crossover from Edwards- Wilkinson 
(EW)22 to Kardar-Parisi-Zhang (KPZ)33 scaling. A rep- 
resentative model in this latter group is the competition 
between ballistic deposition (BD)2i and random deposi- 
tion with surface relaxation (RDSR, or Family model)^, 
respectively with probabilities p and I — p. In all cases, 
the crossover appears for small p. 

For several reasons, the association of those models 
with stochastic growth equations of the Langevin type 
is a problem of central interest. First, it facilitates find- 
ing asymptotic properties which are frequently unclear 
in numerical works on lattice models^^""— . Secondly, 
renormalization study of the growth equation may relate 
unexpected numerical results to crossover or instability 
effects^S. Finally, improvement in atomistic modeling of 
thin film growth may be achieved from the advance on 
the stochastic equation analysis. 

Langevin equations for some of those competitive 
growth models were derived in Refs .*^d^'^° and scaling 
features of KPZ and EW equations were discussed by sev- 
eral authorsi2ii8J^. For RD-RDSR and BD-RDSR mod- 
els, Muraca et ali^ suggested quadratic dependence of 
equation coefficients on the probability p in the crossover 
regime, which was in good agreement with available nu- 
merical data. They argued that the time for the less 
probable process to occur scales as 1/p and that the cor- 
responding increase in local height is also proportional 
to p [Eqs. (1) and (2) of RefJ^]. The resuhs for RD- 
RDSR were confirmed by scaling arguments in Refsi ^'^'d^ . 
However, the claim on the universal quadratic form of 
vanishing coefficients^ is ruled out by the study of a 
restricted solid-on-solid (RSOS)'^l model with deposition 
and erosion, which shows linear p-dependence of the non- 



2 



linear term in the KPZ equation^S, and by the RD-BD 
modeliii^, which shows p^^^ scahng of that term. 

In this work, Langevin-type equations associated with 
various competitive lattice models showing RD-CD and 
EW-KPZ crossovers are derived through a standard van 
Kampen expansion of the Master equation, followed by a 
proper choice of the jump moments. From this approach, 
the form of the equation coefficients is remarkably differ- 
ent from the one proposed by Muraca et aU^: in all cases, 
the random choice of the asymptotically dominant pro- 
cess gives a linear dependence on p for the coefficients 
that vanish as p 0, instead of the quadratic depen- 
dence. However, for RD-CD models, average local slopes 
diverge as p — J> 0, thus the optimal regularization of step 
functions (present in the transition rates of all discrete 
models) have lowest order coefficients that scale as p or 
p^/^, depending on the aggregation mechanism of the CD. 
The combination of those scaling relations give equation 
coefficients of EW and KPZ equations vanishing as p, p^, 
or p'^/^. In all cases, they agree with simulations and 
other scaling approachesiii^. For the RDSR-BD model, 
we show that a non-asymptotic regime with quadratic 
scaling is present, which is associated with the dominant 
effect of subsequent BD events, while the asymptotic lin- 
ear relation is predicted for p much smaller than that of 
previous simulations^'^. 

The rest of this work is organized as follows. In Sec. 
II, three lattice models with the RD-CD crossover are 
defined and the approach to predict amplitudes of rough- 
ness scaling is reviewed. In Sec. Ill, details of the method 
to derive Langevin equations are presented and the equa- 
tion for the RD-RDSR model is obtained. In Sec. IV, 
the equations for the RD-BD model and for the bidis- 
perse ballistic deposition are constructed. In Sec. V, the 
equation for the BD-RDSR model is presented and the 
crossover from quadratic to linear scaling is discussed. 
In Sec. VI we summarize our results and present our 
conclusions. 



II. LATTICE MODELS WITH COMPETITION 
OF CORRELATED AND UNCORRELATED 
DEPOSITION 

In these models, growth begins with a flat d- 
dimensional substrate with cubic symmetry and L ad- 
sorption sites (or columns) in each direction, with a total 
of L'^ sites. Cubic particles of lateral size aj| (parallel 
to the substrate plane) and vertical size a± (parallel to 
the average growth direction) are sequentially released at 
randomly chosen columns above the deposit and fall ver- 
tically towards the substrate. The time interval for depo- 
sition of one layer of atoms \L'^ atoms] is r. Each incident 
particle may irreversibly stick at the top of the column of 
incidence, with probability 1 — p (RD), or move and stick 
following some aggregation rule that takes into account 
the neighboring column heights (mimicking physical pro- 
cesses like diffusion, desorption, or bond formation), with 
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FIG. 1: Microscopic rules of the models (a) RDSR, (b) BD, 
and (c) BBD. The incident particles are shaded squares, the 
deposited particles are empty squares with solid boundaries 
and the squares with dashed boundaries indicate the points of 
first contact with the surface. In RDSR, the incident particle 
may move from the point of first contact, as illustrated by 
the arrows. In the other models, aggregation occurs at that 
point. 



probability p (CD). 

The competitive models where RDSR and BD are the 
correlated components were introduced by Albano and 
co-workersi2iii. The rules of RDSR^^ are illustrated in 
Fig. la: the particle sticks at the top of the column of 
incidence if no neighboring column has a smaller height, 
otherwise it sticks at the top of the column with the 
smallest height among the neighbors (if two or more 
neighbors have the same height, one of them is randomly 
chosen). In BD, which is illustrated in Fig. lb, the in- 
cident particle aggregates at the first site where it finds 
a nearest neighbor occupied site (lateral or below it)^, 
which generates a porous deposit. 

BBD was introduced in Ref.'^^ and is itself defined as a 
competitive model. Taking the d = 1 case for illustration, 
particles of two different sizes incide towards the surface: 
single-site particles (lateral size a\\ , vertical size a_i_) with 
probability 1 — F and double-site particles (dimers with 
lateral size 2a||, vertical size aj_) with probability F . Any 
incident particle permanently sticks at the first position 
where it encounters a previously deposited particle be- 
low it. The aggregation of a dimer leads to pore forma- 
tion when the neighboring columns have different heights, 
similarly to the lateral aggregation of BD. The rules of 
BBD are illustrated in Fig. Ic. 

In these models, the surface configuration is the ar- 
ray of discrete height variables H = {hi\^ where i = 
1, 2, 3, ... , L'^. In all deposits (compact or porous), the 
height variable is the one of the highest particle at that 
substrate site, thus H always refers to the outer surface 
of the deposit. 

The global roughness of the surface is defined as the 
rms fluctuation of the height variables, whose average is 
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h: 



WiL,t) = 




Here the overbars indicate spatial averages over the 
height variables, the sum is taken over all the L'' sub- 
strate sites, and the angular brackets indicate configura- 
tional averages. At short times, RD dominates, thus the 
roughness increases as 

WRD-a^{t/Tf\ (2) 

After a crossover time tc, the CD determines the univer- 
sality class of the process. The roughness follows Family- 
Vicsek (FV) scaling"^^ as 

W{L,t)^AL'^f(^^y (3) 

where a is the roughness exponent, / is a scaling function 
such that / 1 in the regime of roughness saturation 
(t oo) and tx is the characteristic time of crossover to 
saturation, which scales as 

tx ~ BL', (4) 

where z is the dynamic exponent. For t ^ (but t ^ 
tc), the roughness scales as 

W^Ctf^, (5) 

where /3 = a/z is the growth exponent. In this growth 
regime, f{x) ^ in Eq. ([3]). 

The exponents a, /? and z depend on the basic sym- 
metries of the CD, but the amplitudes A, B and C are 
mo del- dependent. For small p, they scale as 



A r- 




B r- 




C r- 





where the convention of crossover exponents {d, y, 7) of 
Albano and co-workers^° was used. FV scaling implies 

y/3-(5 + 7 = 0. (7) 

The scaling approaches of Refs.— explain the values 
of the exponents obtained in simulations of lattice models 
and show that only one exponent [usually y or d; Eq. ([S])] 
is sufficient to characterize the crossover. 

For p ^ 1, most deposited atoms attach to the top of 
the randomly chosen site (RD). Thus, from Eq. the 

height difference of neighboring sites is of order Ahur) ~ 

1/2 

a±{At/T) ' after a time interval At. On the other hand, 
the average time for a correlated event (probability p) to 
occur at a given column is Tc ^ t /p. 

In the case of ballistic-like models (e. g. BD or BBD), 
the frequent lateral aggregation (e. g. at columns 2 and 8 



in Fig. lb) immediately creates correlations between the 
neighboring columns. Thus, the time of crossover from 
random to correlated growth is 

ic ~ - p-^T (8) 

This applies to other BD-like models, as discussed, e.g., 
in Ref.34. 

However, in the case of solid-on-solid (SOS) models, 
such as RDSR, the column height changes by a single 
lattice unit at each time step. A single CD event does not 
cancel the random height fluctuation Ah of neighboring 
columns; instead, it only reduces that fluctuation by one 
lattice unit (e. g. the particle aggregating at column 3 
in Fig. la). The height fluctuation produced by RD will 
be suppressed only when the number of correlated events 
Nc is of order Ahuo / a^. At the crossover time t^ this 
number is Nc = tc/Tc ^ tcp/r, thus 

tc - P'^T. (9) 

In both cases, all time scales of the purely correlated 
system {p = 1) are also changed by the scaling factor 
tc/r, such as the saturation time ix [Eq- Q]- Thus we 
have y — 1 ioT ballistic-like models and y = 2 for SOS 
models with single particle deposition attempts [Eqs. 
(|6])]. This result does not depend on the universality 
class of the CD neither on the substrate dimension. 

The average height difference between neighboring 
columns saturates at 

Ah^ a^itc/rf^ ax_p-y/'^ (10) 

for all models (ballistic-like and SOS). This is the scaling 
factor for global height fluctuations [Eqs. dSj and dH)], 
thus 

5^y/2. (11) 

Combined with Eq. ([7]), it shows that a single exponent 
[y) fully characterizes the crossover. 



III. THE RD-RDSR MODEL 

The equation associated to the RD-RDSR model is 
constructed through a van Kampen expansion of the 
Master equation"^^, as discussed in Refsi^i^^i^. 

The transition rate M^(H';H) from the height config- 
uration YL — {hi] to the configuration H' — {h'^} is 

Win'- H) = i ^ Wk5 {hi,hk +a^)l[6 {h'^,h,), 

k j^k 

(12) 

where the (5-function product represents the condition 
that H and H' differ by the deposition of one only parti- 
cle, and Wk is the rate at which the process hk — s> hk + a± 
occurs. 
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Let K^^'^ and K'y^'' be respectively the first and second 



jump moments of W, computed througt 



H' 



(13) 



and 



Kg\u) = ^(Z^:; - h,)ih'^ - h,)WiU':H). (14) 



H' 



According to a theorem of Kurt a^^'^" , later revisited by 
Fox and Keiser in the context of a macrovariable descrip- 
tion for noisy trajectories^ (see also Ref42,), we expect 
that 



dh 



^-ifW(H) 



gives the macroscopic description of H in the hydrody- 
namic limit. If all the conditions imposed in Ref4^ and 
in Refsi^"— are met, the fluctuations in Eq. ([T5|) must 
obey 



{V^{t)) = 



and 



where we used the identity 



(16) 



(17) 



(18) 



between the first and second jump moments in Eq. (1171) . 
In fact, all higher-order jump moments are proportional 
to k[^\ as can be seen from direct calculation. 

For the RD-RDSR model (and related competitive 
models) in d = 1, the first jump moment [Eq. (|13p ] can 
be cast to the form 



K, 



(1) 



( (0) 1 (1) I (2) ^ 1 n ^' 



(19) 



where each wj'^'' gives the conditions for a particle incid- 
ing at column j to move and stick to one of its neighbours 
or to stick at the incidence column. Those conditions de- 
pend on the local height configuration. 

(k) 

The variables w- are called aggregation rules. For 
the RDSR model, a;'"'' represents the conditions for the 
particle incident at site i to stick at i; wp"*, the conditions 
for the particle to relax to its left site (i — 1); and uji , 
the conditions for the particle to relax to its right site 
(i + 1). The aggregation rules can be written in terms of 
discrete step functions asi^ 



,(0) 



,(1) 



,(2) 



(20) 



where Ol = <d{hj—hk), and Q{x) is the unit step function, 
defined at our convenience to be Q{x) — 1 for x > and 
e{x) = for x < 0. 

In order to pass from the discrete model to its contin- 
uum limit, we assume there exists a continuous function 
'i'{x,t) that interpolates all points hi{t) of the substrate, 
while a|| is kept small but finite. This is possible if we 
can write 



k=l 



dx'^ 



(±a||nj 



(21) 



for some ^'(a;, t). 

We assume further there also exists an analytical rep- 
resentation of the step function Q (see for example 



(15) Refsi^ii^), and we define A^^ 



hk, so that the 



function O(A^) can be expanded in a power series of the 
height differences as 



e(Ai) = 1 



(22) 



where the expansion coefficients have to be chosen 
according to the rules of the lattice model to be 
represented^l"— 

Eq. is inserted in Eq. p3|) . and step functions and 
height differences are expanded according to Eqs. ((22)) 
and (HJ). Retaining terms up to the leading order in a|| 
and a±, and in the limit of small p, we obtain the EW 
equation3^ 



dh 
'dt 



lyV'^h + ri{x,t), 



where 



and 



2a^a| 
V = Aip 

T 



F : 



T 



(23) 



(24) 



(25) 



These coefficients differ from those of Ref J^, which 
gave V = (2a|/T)Aip2 and F = {a^/r) (1 - pf + . 

That work proposes that the height at a given column 
increases by a factor proportional io p {1 — p) after a 
time interval t/p [t/ {1 — p)] characteristic of the RDSR 
(RD) process. This gives the quadratic dependence of 
V on p. However, this hypothesis also leads to a flux F 
depending on p, which is not true. Instead, the model is 
SOS and no deposition attempt is rejected, thus the flux 
is independent of p, as given in Eq. (j25p : one layer of 
atoms of height a±_ is deposited during time t. 

On the other hand, numerical workiS and scaling 
argumenta^ii^ give v ^ p^ for small p, which apparently 
disagrees with Eq. (|M| . As will be explained below, an 
additional p factor is hidden in the coefficient Ai of Eq. 
, which is known to be model-dependent. 
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Since the step function is limited to values and 1, 
the sum of terms in the right-hand side of Eq. is 
expected to be of order 1 or smaller. The step at the 
origin indicates that the first order term AiAj, is finite 
and non-zero. Indeed, when Aj is computed through 
some continuous representation of the step function, as 
in Ref ji^ii^, the first order term is of order 1 for any 
model. 

For pure correlated models, such as RDSR, height dif- 
ferences A-^ are of order 1, thus we sure expect Ai is of 
that order. On the other hand, in RD-RDSR with small 
p, typical neighboring height differences A^. diverge as 
p"^ [Eq. (ITOl) with y = 2]. Thus Ai must vanish as 



(26) 



for a correct regularization of the step function for small 
P- 

Substituting this result in Eq. (1^ , we obtain v ^ p'^. 
Thus, a first factor p in the surface tension coefficient v 
comes from the random choice of correlated events, while 
a second factor p comes from the reduced smoothing ef- 
fect of each correlated event in surface steps of depth 
p"^, as discussed in Sec. II. This interpretation also dif- 
fers from that of Muraca et al^^, that relates the complete 
p^ factor to the random choice of RDSR. 

Analogous arguments apply to the step function ex- 
pansion of RD-CD models in general. From Eq. (ITUl) . 
they give 



A, ^/ 



(27) 



Since the second order term of the expansion in Eq. ([^^ 
must also be of order one, we expect 



Ao 



„25 



' P 



(28) 



Eventually the second-order term in Eq. ([^^ is zero or 
converges to zero as p — >■ 0, which would imply that A2 
is zero or has a higher power in p. 

Another important question is a possible crossover ef- 
fect on the noise term of the stochastic equations. How- 
ever, it can be shown that the noise amplitude is inde- 
pendent of the parameter p. At the end of Sec. IIVI a 
detailed discussion is presented for all models with the 
RD-CD crossover. 



IV. THE RD-BD AND THE BBD MODEL 

For the one-dimensional RD-BD model, the first jump 
moment is 



(1) _ P 



,(4) 



(5), 



+ (1-P)- 



(29) 



with the aggregation rules given by 



CO. 



(3) 



1 



,(4) 



,(5) 



r^ei - ^ (1 



qi \ H+1 

^l+l) "i-l 



(30) 



where (5* — 



J 

,(3) 



— 1 \s the Kronecker delta function. 



Note that a;)"'' + ujf^ + Lof ^ 1 when hi — hi^i or 
hi = /li+i. However, the corresponding jump moment 
expression [Eq. (I^Hl) ] is correct, since oj^'* and oj^-'*'' are 



,(5) 



multiplied by hi^i 

,(3) 



hi and ft-i+i — hi, respectively. Also 
note that a;)"^'' and account for lateral aggregation, 
while refers to aggregation at a local maximum (see 
Fig. 1). 

The same steps of the previous model are then fol- 
lowed: Eq. (|29p is inserted in Eq. (|15|) . with step func- 
tions and height differences expanded according to Eqs. 
((22)) and (|2T|) . Retaining terms up to the leading order in 
a|| and a±, and in the small p limit, we obtain the KPZ 
equation 



dh 
'dt 



F + vV^h 



r]{x,t) 



with 



—P, 

T 



A ^Aip, 



and 



F 



T 



(31) 



(32) 



(33) 



(34) 



A naive inspection of Eqs. ((32|) and (|33p suggests v ^ p 
and X ^ p. However, the scaling of Ai and A2 in Eqs. 
(I?7| and (1^ . with y = 1 for ballistic-like models, gives 



p,\^p 



3/2 



(35) 



which agrees with simulation^ ^ and scaling 
argumentsiiii^. 

Thus, this model also shows that the scaling of the 
equation coefficients on p depend not only on the random 
choice of BD but also on the scaling of height differences. 
Moreover, our analysis show that equation coefficients 
which are noninteger powers or p, such as p^^"^, can be 
predicted by construction of growth equations from the 
microscopic rules, with the noninteger exponent related 
to the step function regularization. 

The coefficient of the surface tension term [Eq. (|32l) ] 
is always positive. Indeed, the lateral aggregation rules 



of BD [u]i ' and w\ ' in Eq. ([5(11) ] reduce local height dif- 
ferences, as illustrated by deposition at columns 2 and 
8 in Fig. lb, which is the role of surface tension. This 
balances the negative contribution to the surface tension 
from aggregation at a local maximum, illustrated by de- 
position at column 5 in Fig. lb. 
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This result differs from what is obtained with the ag- 
gregation rules of Ref.^"^^: v = — (a|ax/T)^ip, which is 

negative for (expected) positive A\. With the aggrega- 
tion rules presented in that work, lateral aggregation does 
not contribute to surface tension, thus aggregation at a 
local maximum renders v negative, as can be inferred 
from the aj^ factor appearing in that formula for v (in 

Ref J^, it follows only from the expression of iJf' ) . 

The preceding discussion shows that it is essential to 
check the consistency of the equation coefficients and the 
geometry of the model when strong approximations (such 
as regularization of step functions) are involved. When 
the surface is rough, pure lateral aggregation at edges 
and steps tends to bring the surface to a plain surface 
state, through both non-conservative (taming of a height 
difference larger than aj_) and conservative mechanisms 
(taming of a step of height aj,). Indeed, we chose aggre- 
gation rules for BD such that the pure lateral aggrega- 
tion gives a positive contribution to the laplacian term, 
as it is expected. That was only possible by allowing 
-f iJ-f^ + iJ'p ^ 1 in Eq. ([5(11) and avoiding products 
between terms such as 1 — 0* , which will be subsequently 
regularized. 

Now we consider the BED model. The first jump mo- 
ment is 



K. 



(1) _ 



(36) 



+ il-p)-, 

T 



with 



OJ. 



(6) ^ 
i 

,(7) _ 



,(8) 



,(9) 



(i/2)r' 

(1/2)0^+1 

(1/2) 

(1/2) (1-^1+1) 



(37) 



where the factors 1/2 correspond to the equal probability 
for the two possible orientations of the dimers on the inci- 
dence site. The KPZ equation [Eq. pip ] is also obtained 
from this rules. In the small p limit, it has coefficients 



—P, 

T 



T 



and 



(l+p). 



(38) 



(39) 



(40) 



The p-dependence of Ai for ballistic-like models [Eq. (p7)) 
with 6 = 1/2] again gives i/ ~ p and A ^ P^^"^, in agree- 
ment with Ref.— , which combined scaling properties of 
the KPZ equation in one dimension and numerical re- 
sults. 

The relation between the first and second jump mo- 
ments (Eq. [T8|) shows that this scaling picture can also 
provide the noise term of the growth equation. We follow 
Ref.i and rewrite Eq. ^ as {7]i{t)Tjj{t')) = D5ijS{t'-t) 

to have D = a±K^^\ where D is the amplitude of the 
noise correlations. For each model, we expand the step 
functions and height differences (Eqs. [22] and [21]) and 
retain terms up to the leading order in a± and a||. This 
gives D a±F for all models, as can be found from in- 
spection of Eqs. ([24]) and ([25]) . ([32]) to ([34]) . ([38]) to pO[) 
and (|44p to ([^5]) , since all time and space derivatives of h 
are finite in the limiting process. That means D = o^/t 
plus terms of order a^a| or a^ay in all RD-CD mod- 
els. For BBD, D = a'i{l+p)/T + ©(a^al), due to the 
particular choice of the time unit for that model. 

These results show that, in the small p limit, there is 
no effect of this parameter on the noise amplitude. Con- 
sequently, all the crossover effects depend on the coeffi- 
cients v and A (in contrast to what is observed in other 
growth models^). 



V. THE RDSR-BD MODEL 

This model was introduced in Ref.^ and involves the 
competition of BD (KPZ class), with probability p, and 
RDSR (EW class), with probability I - p. In Rei^, 
scaling properties were studied numerically, with the co- 
efficient of the nonlinear term scaling as A '~ for 
0.2 < p < 0.5. That quadratic dependence was proposed 
analytically by Muraca et aU^. 

The first jump moment in this case is 



K 



(1) _ 



pIC. 



BD 



RDSR 



where ICf^^^^lCf^ are the first jump moments 



and 



1 r 



RDSR 



r 



,(2) 



(41) 



(42) 



,(3) 



(/i.,_i - hi) + 



(5) 
- cj; a 



(k) 

where w; , k 



(43) 

0, ... 5, are the aggregation rules given 
in Eqs. (|20p and ([30l) . Following the same approach 
of the other models, we obtain the KPZ equation with 
coefficients 



'II, 



2ana_L 



-v4i 1 



2«N 

A = — ^^1 {5-a±Ai)p, 

T 



(44) 



(45) 
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and 



Both RDSR and BD have correlated kinetics which lead 
to finite average values of local slopes, even in the steady 
states. Consequently, the leading coefficients of the step 
functions {Ai, A2) do not vanish in the limit p — )■ 0, in 
contrast with the models with crossover from RD. Thus, 
Eq. (|45t gives X ^ p as p 0, while the other coefficients 
remain nonzero. 

This result disagrees with the quadratic dependence 
observed in simulations of Refji^ and suggested in ana- 
lytical work of Ref J^. In order to understand this dis- 
crepancy, the conditions where the BD component gen- 
erates nonlinearity in the RDSR-BD model have to be 
analyzed. 

In Fig. lb, deposition at columns 2 and 8 shows the 
condition in which lateral aggregation (characteristic of 
BD) leads to excess velocity: deposition occurs at a col- 
umn i which has at least one neighbor with height larger 
by 2a_L or more. This leads to formation of a hole in 
column i. However, for small p, pure RDSR dominates. 
We simulated the one-dimensional RDRS model in lattice 
sizes L — 256 and L — 512 and found that the number of 
columns where excess velocity is possible is P « 0.044. 
This fraction is small because RDSR produces a very 
smooth surface, with a very small number of high steps. 
On the other hand, in pure BD, our simulations show 
that this probabillity is near 1/2. Thus, in the compet- 
itive model with small p, the fraction of columns which 
have lateral growth (i. e. nonlinear growth) is approxi- 
mately Pp. 

On the other hand, the BD model itself creates condi- 
tions for two neighboring sites to have height difference 
2a± or more: if a column i has one larger neighbor j 
{hi < hj), a BD event at j followed by another BD event 
at i leads to lateral aggregation with formation of a hole. 
For instance, this would correspond to the deposition in 
column j — 5, shown in Fig. lb, followed by deposition 
in column i = 4 (not shown) . In the pure RDSR surface, 
the fraction of columns with at least one higher neighbor 
is Q w 0.44, also obtained from simulation. Thus, the 
probability of this column having nonlinear growth due 
to those subsequent BD events is approximately Qp^. 

From the point of view of Eq. (P5|) , pure RDSR corre- 
sponds to Ai P in the regularization of step functions, 
while a BD event corresponds to Ai ^ Qp. 

For small enough p, we certainly have Pp > Qp^ , thus 
the crossover EW-KPZ is dominated by BD events tak- 
ing place on a nearly pure RDSR surface. This occurs for 
p < P/Q ~ 0.10. In this regime, the linear dependence 
of A on p [Eq. (j45p ] is expected. However, simulation re- 
sults of Ref.^'* are for p > 0.2. In the lower limit p — 0.2, 
we have Pp « 0.0088 and Qp"^ « 0.0176. This means 
that Qp^ is twice as large as Pp, and their difference 
is enhanced for larger p. Consequently, the simulated 
range of p favors nonlinearities arising from two subse- 
quent BD events at neighboring columns, which explains 



the observed quadratic dependence of the coefficient A on 
P- 

The arguments of Ref for the X ^ p^ behavior were 
also based on the association of the p^ factor to the ran- 
dom choice of BD. However, it is also a double counting of 
the factor p, which is reasonable for the simulated range 
of p but fails for very small p. 

The linear scaling of the coefficient A on p was also 
found in the RSOS model of deposition and erosion of 
Rcf.'''^, both in simulations and in the derivation of the 
associated KPZ equation. The constraint on the neigh- 
boring height difference of the RSOS model leads to re- 
jection of deposition and erosion attempts, which is the 
mechanism to generate nonlinear growth. That rejection 
occurs with a probability much larger than the probabil- 
ity P for the RDSR, thus the linear dependence on p was 
easily observed in simulations^^. 



VI. CONCLUSION 

Langevin equations associated with various competi- 
tive lattice models were derived. The approach is based 
on a van Kampen expansion of the Master equation, but 
the correct assessment of how do characteristic times and 
lengths scale with the competing parameter plays a cen- 
tral role if we are to find the true dependence of the 
equation coefficients in the crossover regimes. Moreover, 
it is essential to choose representations of aggregation 
rules (using e. g. step and delta functions) that lead to 
physically reasonable equation coefficients, as RDSR and 
BD models illustrate. 

We considered a series of models with crossover from 
random deposition to correlated growth (RD-CD), with 
probability p for the latter, and a model with EW-KPZ 
crossover, with probability p for the KPZ component. 
All coefficients that vanish as p — show a linear p de- 
pendence arising from the random choice of aggregation 
rules. However, in the RD-CD case, neighboring height 
differences diverge in that limit, which leads to the p- 
scaling of the parameters of the optimal regularization of 
step functions. Thus, the coefficients depending on those 
parameters show scaling as p, p^, and p^^^, in all cases 
in agreement with simulation results and other scaling 
approaches. For the model with EW-KPZ crossover, the 
quadratic dependence of the nonlinear term coefficient, 
observed in simulations, is explained as a crossover be- 
havior due to particular model features, while linear p 
dependence is expected for very small p. Although the 
scaling properties derived here are similar to previous 
works on those model a^^i^^'^^ , the interpretation is very 
different and the applicability of the method is broader, 
for instance being extendable to higher dimensions. 
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